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ABSTRACT 

The assembly of galaxies can be described by the distribution of their star formation as a function of cosmic time. 

Thanks to the WFC3 grism on the Hubble Space Telescope ( HST) it is now possible to measure this beyond the 
local Universe. Here we present the spatial distribution of Ha emission for a sample of 54 strongly star-forming 
galaxies at z ~ 1 in the 3D-HST Treasury survey. By stacking the Ha emission, we find that star formation occurred 
in approximately exponential distributions at z ~ 1, with a median Sersic index of n = 1.0 ± 0.2. The stacks are 
elongated with median axis ratios of b/a = 0.58 zb 0.09 in Ha consistent with (possibly thick) disks at random 
orientation angles. Keck spectra obtained for a subset of eight of the galaxies show clear evidence for rotation, with 
inclination corrected velocities of 90-330 km s -1 . The most straightforward interpretation of our results is that star 
formation in strongly star-forming galaxies at z ~ 1 generally occurred in disks. The disks appear to be “scaled-up” 
versions of nearby spiral galaxies: they have EW(Ha) ~ 100 A out to the solar orbit and they have star formation 
surface densities above the threshold for driving galactic scale winds. 

Key words: galaxies: evolution - galaxies: formation - galaxies: high-redshift - galaxies: kinematics 
and dynamics - galaxies: star formation - galaxies: structure 
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1. INTRODUCTION 

Galaxy formation is a complex process involving starbursts, 
mergers, and strong gas flows (e.g., Hopkins et al. 2006; Brooks 
et al. 2009; Dekel et al. 2009). Furthermore, stellar migration 

and secular processes can change the structure of galaxies at late 
times (e.g., Roskar et al. 2008; Grand et al. 2012). Therefore, 
even if we could perfectly locate and age date every star in the 
Milky Way, we still could not say where and with what structural 

and kinematic properties those stars formed. The only way to 
establish where a galaxy’s stars formed, and hence how it was 

assembled, is to map the star formation while those stars were 
forming. 

Obtaining Ha and stellar continuum maps of z ^ 1 galaxies, 
with the ~ 1 kpc resolution necessary to put constraints on the 

spatial distribution of star formation, is challenging and has so 
far only been possible using a combination of adaptive optics 
and integral field units on 8-10 m class telescopes. These studies 
paint a complex picture: they find that star-forming galaxies at 
z ~ 2 are a mix of “puffy” and often clumpy rotating disks, 
mergers, and more compact dispersion-dominated objects (e.g., 
Genzel et al. 2008; Shapiro et al. 2008; Cresci et al. 2009; Forster 

Schreiber et al. 2009, 2011; Law et al. 2009; Jones et al. 2010; 
Mancini et al. 2011). 

Some studies claim that there are trends in the structural 
properties of Ha as a function of redshift, implying that the way 
galaxies assemble their stars varies fundamentally as a function 


of cosmic time. In particular, Epinat et al. (2009) and Kassin 
et al. (2012) suggest that galaxies become cooler and more 
rotation-dominated with time, z ~ 1-0 being the epoch of “disk 
settling.” This is interesting because most of the stars in the 
disks of galaxies like the Milky Way were formed in this epoch. 

It is now possible to obtain high spatial resolution (<1 kpc) 
information on Ha emission at z ~ 1 with a high Strehl ratio, 
owing to the near-IR slitless spectroscopic capabilities provided 
by the WFC3 camera on the Hubble Space Telescope (HST). 
In Nelson et al. (2012), we used data taken as part of the 
3D-HST survey to build on previous ground-based studies by 
mapping the Ha and stellar continuum with high resolution for 
a sample of 57 galaxies at z ~ 1 and showed that star formation 
broadly follows the rest- frame optical light, but is slightly more 
extended. Here we stack the Ha maps of these galaxies to 
construct high signal-to-noise ratio (S/N) radial profiles and 
measure structural parameters of stellar continuum emission 
and star formation. These unique, high spatial resolution data 
from 3D-HST are combined with kinematics from the Near 
Infrared Spectrometer (NIRSPEC) on the W. M. Keck telescope 
(McLean et al. 1998). We assume H 0 = 70 km s _1 Mpc -1 , 
Q m = 0.3, and Q A = 0.7. 

2. 3D-HST AND SAMPLE SELECTION 

The 3D-HST survey, a 248-orbit Treasury program on the 
HST , supplies the two-dimensional emission line maps needed 


1 


The Astrophysical Journal Letters, 763:L16 (6pp), 2013 January 20 


Nelson et al. 


>% 

o 


+ 

> 

D 

q: 

LL 

c d 

CD 

O 


3 


2 


1 


0 


-1 



9 10 11 

log(M) [MJ 


Figure 1. Location of selected galaxies in the SFR-mass plane. Six sources 
without IR photometry are not shown in the figure. The histogram shows the 
position of the selected galaxies relative to the rest of the galaxies by dividing 
a power-law fit out of the distribution. The purple line shows the fit, the light 
purple region delineates ±lcr. The selected galaxies appear to lie largely on the 
upper half of the “main sequence” in SFR(IR+UV). 

(A color version of this figure is available in the online journal.) 


to measure the spatial distribution of star formation. The 
WFC3 G141 grism provides spatially resolved spectra of all 
sources in the field. The wavelength range of the G141 grism, 
1.15 /rm < A. < 1.65 /rm, covers the Ha? emission line for 
galaxies in the redshift range 0.7 < z < 1.5. The survey also 
provides broadband near-infrared imaging in the F140W filter, 
which samples the rest-frame R band for z ~ 1. Using the 
same camera under the same conditions to map both Ha and 
continuum emission allows us to compare their distributions 
directly, a crucial advantage of this strategy. The maps have 
a spatial sampling of ~0.5 kpc at z ~ 1 (0?06 pixels). The 
data presented in this Letter are based on the ~70 pointings, 
roughly half of the full data set, that were obtained prior to 201 1 
June (see Nelson et al. 2012). The survey, data reduction, and 
determination of derived quantities (e.g., redshifts and masses) 
are described in Brammer et al. (2012) and van Dokkum et al. 
(2011). 

The 3D-HST grism spectra have high spatial resolution and 
low spectral resolution ( R ~ 130), meaning that emission line 
structure reflects almost exclusively spatial structure (morphol- 
ogy) in contrast to data with high spectral resolution where 
structure reflects velocity (rotation or dispersion). Emission line 
maps of galaxies are made by subtracting the continuum emis- 
sion from the two-dimensional spectrum. The details of the 
procedure are described in Nelson et al. (2012), Lundgren et al. 
(2012), and K. B. Schmidt et al. (2013, submitted). 

We selected galaxies with redshifts 0.8 < z < 1.3, total 
F140W AB magnitude <21.9, and rest-frame EW(Ha?) > 
100 A. We also limited the sample to galaxies with effectively 
no contaminating flux from the spectra of other nearby objects. 
The position of the sample in the 0.8 < z < 1.3 star formation 
rate (SFR)-stellar mass plane is shown in Figure 1 . SFRs were 
derived from the IR + UV luminosities as described in Whitaker 
et al. (2012). The selected galaxies lie predominantly on the 
upper half of the “star-forming main sequence” (as shown in 


the histogram of Figure 1; e.g., Noeske et al. 2007; Daddi et al. 
2007; Whitaker et al. 2012). These galaxies were selected in 
Ha not IR; while this selection may bias the sample, it may 
possibly also have the advantage that it adds less uncertainty to 
interpreting the spatial distribution of star formation due to dust 
extinction. 

3. THE SPATIAL DISTRIBUTION OF STAR FORMATION 
3.1. Stacked Ha and Stellar Continuum Emission 

We stack the high spatial resolution Ha maps from 3D-HST 
to create average Ha maps — increasing the S /N and providing 
for a reliable measurement of the structural properties of Ha to 
large radii for this sample at z ~ 1 . We divide the sample in three 
bins based on their Ho? size 10 in order to homogenize the sample 
and investigate trends with size. * 11 The median stellar masses of 
the stacks are log(M*) = 10.0, 10.2, 10.2. The galaxies are 
normalized by their F140W flux and centered according to 
their luminosity- weighted F140W image centers. We masked 
the [Sn] AA6716, 6731 A emission in the fit. At our spectral 
resolution, the FWHM of a line is ~100 A, meaning that Ha 
A.6563 A and [Nil] A6583 A are unresolved but Ha and [Sn] 
are separated by ~3 resolution elements. We rotated the images 
based on their GALFIT-derived F140W position angles to align 
them along their major axes, summed them, and divided the 
resulting stack by the summed masks. Finally, in order to correct 
for the effects of the point-spread function (PSF), the stacks were 
deconvolved by adding the GALFIT-derived PSF-deconvolved 
model to the residuals of the fit, a method developed in Szomoru 
et al. (2010). 

Figures 2 and 3 show the stacked, PSF-corrected Ha emis- 
sion. The S /N in the stacks is high, and the emission is clearly 
spatially extended. Furthermore, both the rest- frame R band and 
Ha maps are clearly elongated. We quantified the structural 
properties of the stacks with GALFIT (Peng et al. 2002) using 
empirical PSFs. For the continuum and Ha?, we find weighted 
mean axis ratios of b/a(R ) = 0.55 zb 0.07 and b/a(Ha) = 
0.58 =b 0.09, respectively, consistent with expectations for ran- 
domly oriented disks with thickness c/a = 0.26, 0.33. This 
thickness could be attributed to disks being intrinsically thick at 
z ~ 1, heterogeneous properties of the galaxies being stacked, 
the effect of bulges, and a selection biased toward face-on galax- 
ies, or noise. 

3.2. Radial Profiles are Nearly Exponential 

Figure 3 shows the radial profiles derived from the stacks. 
Radial profiles of the PSF-corrected stacks were created by 
summing the flux in elliptical radii. A 30% correction was 
applied to account for [N ii] (see Section 4). All uncertainties 
were derived by bootstrap resampling the stacks. 

The Ha? emission ranges from exponential to much less 
centrally concentrated than exponential. The stellar continuum 
emission in all the stacks is nearly exponential. 

The radial equivalent width profiles (right panels of Figure 3) 
are flat to within a factor of two and EW(Ha?) > 100 A 
for all stacks within 1.5 r e . This suggests that the high global 
EW(Ho?)s seen in these galaxies are not driven primarily by 
central starbursts but by strong star formation at all radii. 


10 The sizes are measured using growth curves as described in Nelson et al. 
( 2012 ). 

11 The results are similar when the bins are defined according to the F140W 
sizes. 
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Figure 2. Stacks of the high spatial resolution, PSF-corrected maps of Ha emission (top). Stacking was done based on Ha size; the number of galaxies included 
in each size bin is listed (Attack) as is the major axis effective radius (in kpc) of each stack measured by galfit ( r e ). The stacks have a weighted mean Sersic index 
n{ Ha) ~ 1 and axis ratio of b/a{ Ha) = 0.58 ± 0.09, consistent with disks at random orientation angles. Bottom panels show corresponding radial profiles with dark 
purple — small, medium purple — mid-sized, and light purple — large (Figure 3, Section 3.2). The stacked Ha emission always has Sersic index n < 1. 

(A color version of this figure is available in the online journal.) 







Figure 3. Stacks of Ha (left) and rest-frame R-band (middle) emission have nearly exponential (or shallower) radial profiles. The EW(Ha) profile (discussed in 
Section 3.2) is shown in the right panel. The bottom panels show the profiles normalized by their effective radius. The horizontal line is the SFR surface density 
criterion for driving large-scale outflows (Heckman 2002). As in Figure 2, the dark, medium, and light colors correspond to the small, medium, and large stacks, 
respectively, as shown at the top. The spatial distribution of Ha emission is exponential or shallower. 

(A color version of this figure is available in the online journal.) 
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Figure 4. Left panel shows the radial color profiles for each of the stacks with black for the small stack, dark gray for the medium stack, and light gray for the large 
stack. The right panel shows the implied dust corrected radial profiles of star formation surface density with the colors as in Figures 2 and 3. Error bar on the right 
panel denotes a typical uncertainty. Dust- corrected radial profiles of star formation remain nearly exponential. 

(A color version of this figure is available in the online journal.) 


The derived Sersic indices (Figure 2) of the Ha and stellar 
continuum profiles are n(Ha) = 1 .2 ± 0.6, 0.9 ±0.1, n(Ha) < 
1 andft(F140W) = 1.8 ± 0.4, 1.4 ±0.1, 1.1 ± 0.2 with weighted 
means of n(Ha) = 1.0 ± 0.4 and rc(F140W) = 1.4 ±0.2. As 
shown in van Dokkum et al. (2010), structural parameters mea- 
sured from stacks are close to the mean values for the individual 
galaxies going into the stacks. The upper limit of n(Ha ) < 1 
for the largest Ha stack reflects the fact that the derived n(Ha) 
depends somewhat on the details of the fit (treatment of sky, 
fitting region) but always n(Ha ) < 1. 

We find that the stacked spatial distribution of Ha for galaxies 
in this sample is exponential or shallower. All stacks have n < 2, 
implying a bulge fraction of less than 20% (van Dokkum et al. 
1998). The radial Ha profiles (Figure 3, left), seem to show a 
trend toward lower n(Ha) with increasing size. However, this 
trend is not statistically significant given the errors in the derived 
Sersic indices. Additionally, note that the median galaxy with 
EW(Ha) > 100 A has a high enough star formation surface 
density to drive winds out to ~1 r e (Heckman 2002). 

3.3. Effects of Dust 

A major uncertainty in the interpretation of the radial Ha 
profile is the effect of differential dust extinction: if some parts 
of the galaxies are more obscured than others (see e.g., Wuyts 
et al. 2012), the derived radial profile of Ha would not reflect 
the radial profile of star formation. To assess the importance of 
this effect, we estimate the extinction as a function of radius. 
We determine the extinction from the radial color profile. 

The color profile is determined from the combination of the 
F140W stacks with ACS F814W stacks, which can be converted 
into rest- frame U—V color profiles. As shown in the left panel 
of Figure 4, for most of the radial extent of the stacks the color 
is fairly constant, but it becomes redder inside a radius of 2 kpc 
(see also Szomoru et al. 2012). To estimate roughly how color 
translates into missed star formation, we assume that the star 
formation that is not captured by Ha will be captured by the IR 
(Kennicutt et al. 2009). Using photometry from R. E. Skelton 
et al. (2013, in preparation), in addition to rest- frame U—V 
colors and IR-based SFRs from the NEWFIRM Medium Band 
Survey (Whitaker et al. 2012), empirical relations were derived 
for the translation of rest-frame U—V colors into star formation 


corrections: 

log(SFR(IR+Ha)) - log(SFR(Ha)) - 0.3 x (U - V). (1) 

The dust-corrected radial profiles are shown in the right panel 
of Figure 4. The smallest galaxies appear to have star formation 
that is marginally steeper than exponential. The large galaxies 
have the largest radial gradients but their implied star formation 
is still less steep than exponential. Although this analysis assigns 
the entire observed color gradient to extinction, the centers of 
the radial color profiles could be red because of dust or age. 
Importantly, the Sersic indices of the star formation in these 
stacks based on the implied dust correction remain close to 
one: n = 1.4, 1.1, 0.6 for each of the stacks, weighted mean = 
1.0 ± 0.4. So, even when accounting for dust, the averaged 
star formation in these galaxies has a nearly exponential spatial 
distribution. 

4. KINEMATICS 

The flattening of the stacks and the exponential Ha profiles 
suggest that the star formation occurs in rotating disks. To test 
this, we measured kinematics for a subset of this sample using 
NIRSPEC on the Keck II telescope on 2012 April 9-10. The 
sample comprises eight galaxies chosen from the small (1), 
middle (5), and large (2) stack, which have sizes, masses, Sersic 
indices, SFRs, and axis ratios representative of the sample as 
a whole. 

We used the low dispersion mode of NIRSPEC with O'/ 5 
seeing and a slit width of 0/7, giving a spectral resolution of 
a ~ 80 km s -1 in the J band (compared to ~500 km s -1 in the 
grism spectra). The slit was aligned along the major axis of each 
galaxy and observations were conducted in a series of four 900 s 
exposures, dithering along the slit. The data reduction followed 
standard procedures for long-slit spectroscopy (see, e.g., van 
Dokkum et al. 2004). We extracted kinematic information from 
the two-dimensional spectra by fitting a Gaussian to the Ha and 
[Nil] emission simultaneously at each spatial position along 
the slit. The median [Nil] /Ha is 0.3. Assuming the velocity 
shear in the two-dimensional spectra is due to rotation, we take 
the rotational velocity to be the velocity difference between the 
geometrical center of the galaxy and the maximum velocity. 
We correct the rotation velocities for inclination angle using 
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Figure 5. Shown here are the rotation curves (middle row) derived from the NIRSPEC spectra (insets in corners), one example from each stack with the implied 
rotation velocities (p CO rr — corrected and u unC orr — uncorrected for inclination) and velocity dispersion (cx) listed in km s _1 . The top row shows the corresponding 
rest-frame R-band images and the bottom row shows a false color image with the Ha emission in red and stellar continuum in blue. In these rows, ellipses mark the 
R-band and Ha effective radii (r e (R), r e (Ha)), respectively, and gray arrows have a scaled length corresponding to the projected specific angular momentum. 

(A color version of this figure is available in the online journal.) 


the 3D-HST axis ratios. We calculate one-dimensional velocity 
dispersions by fitting a Gaussian to the Ha emission in the 
central row of each spectrum. We calculate the instrumental 
a to be ~80 km s -1 by fitting a Gaussian to the sky lines. 
The intrinsic velocity dispersion is calculated by subtracting 
the instrumental dispersion from the measured dispersion in 
quadrature. 

We find that all galaxies show velocity shear in their 
two-dimensional spectra, with derived rotation velocities of 
90-235 km s -1 uncorrected and 110-330 km s -1 corrected for 
inclination. 12 Figure 5 shows the rotation curves for one exam- 
ple galaxy for each stack. These spectra are of the best quality 
but appear to be representative of the eight. The velocity profile 
of the smallest galaxy in Figure 3 does not show evidence for 
a turnover, which means the measured maximum velocity is a 
lower limit on the rotation velocity at large radii. We conclude 
that the kinematics are consistent with the disk interpretation of 
the structural properties of the Ha emission. Although different 
classification methods make it difficult to perform a quantita- 
tive comparison, our results are qualitatively consistent with the 
finding that a large fraction of star-forming galaxies at z ~ 1 


12 Measured velocity dispersions are also relatively high, but are difficult to 
interpret given the low spatial resolution of the Keck spectra. 


appear to be rotating, (e.g., Wisnioski et al. 2011; Epinat et al. 
2012; Swinbank et al. 2012). 

5. DISCUSSION 

The central result of this Letter is that the radial distribution of 
stacked Ha emission in z ~ 1 galaxies is close to exponential out 
to ~10 kpc. Combined with the axis ratios of the stacks and the 
kinematics of a subset of the sample, the most straightforward 
interpretation is that star formation seems typically to occur in 
disks at z ~ 1 at least for galaxies with high EW(Fla'). The 
factor of -10 increase in the SFRs of galaxies from z = 0 
to z = 1 (e.g., Damen et al. 2009; Fumagalli et al. 2012) is 
apparently driven by increased star formation activity in disks 
rather than a much greater prevalence of merger-driven central 
starbursts — consistent with other studies (e.g., Rodighiero et al. 
2011; Wuyts et al. 2011,2012). 

This result raises a number of questions. First, it is not clear 
how these galaxies are related to z ~ 0 galaxies. An average 
galaxy in this sample has a distribution of stars with a Sersic 
index of n = 1.4. If this average galaxy forms stars with n = 1.0, 
it will have a lower Sersic index at later times. Either we are 
witnessing the build-up of only the latest of late-type galaxies 
(these galaxies substantially change where their star formation is 
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occurring between z ~ 1 and z ~ 0) or their stars, once formed, 
must migrate into a different configuration. In other words, if 
these are the ancestors of typical z ~ 0 spiral galaxies, their 
bulges need to be built in some way besides the star formation 
we are seeing. This could be accomplished by mergers or secular 
evolution (e.g., Bournaud et al. 2009; Ellison et al. 201 1). These 
observations could be understood in the context of the following 
picture: at z > 1, disks are gas-rich and turbulent. As redshift 
decreases, the gas fraction decreases, and a bar instability grows, 
increasing the central concentration and central vertical velocity 
dispersion, (e.g., van den Bosch 2001; DeBuhr et al. 2012; 
Forbes et al. 2012). 

We also note the somewhat surprising presence of a popula- 
tion of large, rapidly rotating disks with relatively low stellar 
masses at z ~ 1. These galaxies have mean r e (Ha) = 7.3 kpc, 
v = 240 km s -1 , a = 89 km s -1 , and M* = 2.2 x 10 10 M 0 , 
meaning they are disks 8 Gyr ago, larger in size than the Milky 
Way (Drimmel & Spergel 2001), and thicker (Freeman 1987), 
with ~l/3 of the stellar mass (McMillan 2011) and similar or 
higher maximum rotation velocities (Bovy et al. 2009). Both 
what these galaxies become in the local universe and how such 
extended star-forming disks were made in an epoch of high 
disk turbulence are open questions. The small disks on the other 
end of the size distribution are also of interest. With median 
r e ( Ha) = 1.9 kpc, r e (R) = 2.1 kpc, M* = 8.9 x 10 9 M 0 , and 
stacked n( Ha) = 1.2, they more closely resemble compact ver- 
sions of disky star formers than merger-driven star formation in 
spheroids (see also van der Wei et al. 2011). 

There are several caveats. First, our sample is not complete 
in mass or in SFR(UV+IR). In future papers (using the full 
3D-HST survey) we will study the distribution of star formation 
as a function of these parameters. Second, the distribution of 
star formation in the stacks is not necessarily representative of 
individual galaxies. In individual galaxies the star formation 
is clumpy (e.g., Genzel et al. 2008, 2011; Forster Schreiber 
et al. 2009; Nelson et al. 2012; Figure 1) and it is difficult 
to quantify the structure. Stacking these clumpy objects could 
be hazardous, but it may also provide more insight than can 
be gleaned from individual galaxies: assuming that the clumps 
are transient and “light up” a part of the underlying gas disk 
for a short time (as in, e.g., Wuyts et al. 2012), our stacking 
technique effectively produces a time-averaged map of the star 
formation in the galaxies. Finally, dust attenuation remains a key 
uncertainty, both in the selection and in the interpretation of the 
profiles. This can be addressed with maps of the IR emission at 
the full AFMA resolution, or by measuring (spatially resolved) 
B aimer decrements. 

We thank the referee for a thorough report which improved 
the Fetter. Support from grant HST GO- 12 177 is gratefully 
acknowledged. 
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